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In this paper we derive an extra class of non-Markovian master equations where the system 
state is written as a sum of auxiliary matrixes whose evolution involve Lindblad contributions 
with local coupling between all of them, resembling the structure of a classical rate equation. The 
system dynamics may develops strong non-local effects such as the dependence of the stationary 
properties with the system initialization. These equations are derived from alternative microscopic 
interactions, such as complex environments described in a generalized Born-Markov approximation 
and tripartite system-environment interactions, where extra unobserved degrees of freedom mediates 
the entanglement between the system and a Markovian reservoir. Conditions that guarantees the 
completely positive condition of the solution map are found. Quantum stochastic processes that 
recover the system dynamics in average are formulated. We exemplify our results by analyzing the 
dynamical action of non-trivial structured dephasing and depolarizing reservoirs over a single qubit. 
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I. INTRODUCTION 

The description of open quantum systems in terms of 
local in time evolutions is based in a weak coupling and 
Markovian approximations [l], Q ■ When these approxi- 
mations are valid, the dynamics can be written as a Lind- 
blad equation [U, 0, 0, 13| ■ The evolution of the density 
matrix Ps{t) of the system of interest reads 

= ^[H,ff,ps{t)]-{D,psit)U + F[psit)], (1) 

where i?e// is an effective Hamiltonian, {■••}+ denotes 
an anticonmutation operation, and 

D = IY1 ^o.-rV^^^> = E ^o,yV^ ' ^7 • (2) 



Here, the sum indexes run from one to (dimWg)^, where 
dim 7^5 is the system Hilbert space dimension. The set 
{Va} corresponds to a system operator base, and Ua^y de- 
notes a semipositive Hermitian matrix that characterize 
the dissipative time scales of the system. 

Outside the weak coupling and Markovian approxima- 
tions, it is not possible to establish a general formalism 
for dealing with non-Markovian system-environment in- 
teractions [i, i, S S, S [3, [H El- Nevertheless, there 
exist an increasing interest in describing open quantum 
system dynamics in terms of non-Markovian Lindblad 

equations [H, M M,MM,M, M M 

Here, the density matrix pg{t) of the system evolves as 



dt 



[Heff,Ps{t)]+ dTK{t-T)C[ps{T)], (3) 



'present address 



where £.[•] = —{!?,•}+ -I- is a standard Lindblad 
superoperator. The memory kernel K{t) is a function 
that may introduces strong non-Markovian effects in the 
system decay dynamics. 

The study and characterization of this kind of dynam- 
ics is twofold: on one hand, there is a general funda- 
mental interest in the theory of open quantum systems 
to extend the methods and concepts well developed for 
Markovian dynamics to the non-Markov case. On the 
other hand there are many new physical situations in 
which the Markov assumption, usually used, is not ful- 
fill and then non-Markovian dynamics has to be intro- 
duced. Remarkable examples are sin gle fluorescent sys - 
tems hosted in complex environments [24 [25l. \2T\. \2§^ . 
superconducting qubits [29l [30| and band gap materials 
[3l|,|3l|. 

Most of the recent analysis on non-Markovian Lindblad 
evolutions [H, [fl [H, M, IH, [H were focus on the 
possibility of obtaining non-physical solution for Pplt) 
from Eq. This problem was clarified in Refs. [14l.[l5j. 
where mathematical constraints on the kernel K(t) that 
guarantees the completely positive condition [1, H, 13 
the solution i nap ps{0) Psi^) were found. Further- 
more, in Ref. [l5| the completely positive condition was 
associated with the possibility of finding a stochastic rep- 
resentation of the system dynamics. 

There also exist different analysis that associate evolu- 
tions like Eg. with microscopic system environment 
interactions [11 H El, 113] • In Ref. [U the micro- 
scopic Hamiltonian involves extra stationary unobserved 
degrees of freedom that modulate the dissipative coupling 
between the system of interest an a Markovian environ- 
ment. This kind of interaction lead to a Lindblad equa- 
tion characterized by a random rate. A similar situation 
was found in Ref. [22| by considering a complex envi- 
ronment whose action can be described in a generalized 
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Born-Markov approximation (GBMA). This approach re- 
hes in the possibihty of spUtting the environment in a 
"direct sum" of sub-reservoirs, each one being able to 
induce by itself a Markovian system evolution. When 
the system-environment interaction does not couples the 
different subspaces associated to each sub-reservoir, the 
system dynamics can also be written as a Lindblad equa- 
tion with a random dissipative rate. After performing 
the average over the random rate, the system dynamics 
can be written as a non-local evolution with a structure 
similar to Eq. ([3]). Besides its theoretical interest, the 
GBMA was found to be an useful tool for describing and 
modeling specific physical situations, such as the fluores- 
cence signal scattered by individual nanoscopic systems 
host in condensed phase environments [28l |. 

The aim of the present work is to go beyond previ- 

ous results [H, M M, M [13, H El, H EH SI HI, 

and present an alternative kind of evolution that in- 
duces strong non-local effects, providing in this way an 
extra framework for studying and characterizing non- 
Markovian open quantum system dynamics. In the 
present approach, the system density matrix can be writ- 
ten as 

R 

where the unnormalized states Puit) have associated an 
effective Hamiltonian H'^J^^ , and their full evolution is 
defined by 

- ^[H^d^ ,~PR{t)] - {Dn,~PR{t)]++Fn[Mi^] 

- E {DBJB„~PR{t)}++ E FRR,[p^,{t)l 

r! r! 

R'^R R'^R 

(5) 

subject to the initial conditions 

p^(0) = PrPs{0)- (6) 

The positive weights Pr satisfy Ei? Pr = 1- the other 
hand, the diagonal superoperator contributions are de- 
fined by 

a, 7 a, 7 

while the non-diagonal contributions reads 

Dr'r - ^E«M^- ^««'H = E«fi«'^" • ^7- 

a, 7 a, 7 

(8) 

By convenience, we have introduced different notations 
for the diagonal and non-diagonal terms. As in standard 
Lindblad equations, Eq. ([1]), the matrixes a'^ and a'^j^ 
characterize the dissipative rate constants. The structure 
of the non-diagonal terms in Eq. ^ resemble a classical 
rate equation [33j. Therefore, we name this kind of evo- 
lution as a Lindblad rate equation. 



Our main objective is to characterize this kind of equa- 
tions by finding different microscopic interactions that 
leads to this structure. Furthermore, we find the con- 
ditions that guarantees that the solution map ^5(0) — > 
Ps{t) is a completely positive one. 

While the evolution of pg{t) can be written as a non- 
local evolution [sec Eq. ([CT]) ]. the structure Eq. (O leads 
to a kind of non-Markovian effects where the station- 
ary properties may depend on the system initialization. 
In order to understand this unusual characteristic, as in 
Ref. [1^, [l^l , we also explore the possibility of finding a 
stochastic representation of the system dynamics. 

We remark that specific evolutions like Eq. ([5]) were 
derived previously in the literature in the context of dif- 
ferent approaches [lOj IHi ■ The relation between those 
results is also clarified in the present contribution. 

The paper is organized as follows. Is Sec. II we derive 
the Lindblad rate equations from a GBMA by considering 
interactions Hamiltonians that has contribution terms 
between the subspaces associated to each sub-reservoir. 
An alternative derivation in terms of tripartite interac- 
tions allows to find the conditions under which the dy- 
namic is completely positive. A third derivation is given 
in terms of quantum stochastic processes. In Sec. Ill we 
characterize the resulting non-Markovian master equa- 
tion. By analyzing some simple non-trivial examples that 
admits a stochastic reformulation, we explain some non- 
standard general properties of the non-Markovian dy- 
namics. In Sec. IV we give the conclusions. 

II. MICROSCOPIC DERIVATION 

In this section we present three alternative situations 
where the system dynamics is described by a Lindblad 
rate equation. 

A. Generalized Born-Markov approximation 

The GBMA applies to complex environments whose 
action can be well described in terms of a direct sum 
of Markovian sub-reservoirs [13]. This hypothesis im- 
plies that the total system-environment density matrix, 
in contrast with the standard separable form [l|, [§] , as- 
sumes a classical correlated structure 4j (see Eq. (6) in 
Ref. (2^). In our previous analysis, we have assumed 
a system-environment interaction Hamiltonian that does 
not have matrix elements between the subspaces associ- 
ated to each sub-reservoir. Therefore it assumes a direct 
sum structure (see Eq. (5) in Ref. [13]). By raising up 
this condition, i.e., by taking in account arbitrary inter- 
action Hamiltonians without a direct sum structure, it 
is possible to demonstrate that the GBMA leads to a 
Lindblad rate equation, Eq. ([5]). 

As in the standard Born-Markov approximation, the 
derivation of the system evolution can be formalized in 
terms of projector techniques [ll|. In fact, in Ref. [12] 
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Breuer and collaborators introduced a "correlated pro- 
jector technique" intended to describe situations where 
the total system-environment density matrix does not as- 
sume an uncorrelated structure. Therefore, the system 
dynamics can be alternatively derived in the context of 
this equivalent approach. The main advantage of this 
technique is that it provides a rigorous procedure for ob- 
taining the dynamics to any desired order in the system- 
environment interaction strength [O, [l2|. Here, we as- 
sume that the system is weakly coupled to the environ- 
ment. Therefore, we work out the system evolution up 
to second order in the interaction strength. 

We start by considering a full microscopic Hamilto- 
nian description of the interaction of a system S with its 
environment B 



Ht — Hq + Hb + Hi 



(9) 



The contributions Hs and Hb correspond to the system 
and bath Hamiltonians respectively. The term Hi de- 
scribes their mutual interaction. 

The system density matrix follows after tracing out 
the environment degrees of freedom, Pg{t) = TrB{pj.(i)}, 
where the total density matrix prp{t) evolves as 

^P^ = ^[HT.Pr{t)]^CT[PT{t)]. (10) 
Now, we introduce the projector V defined by 



TrsiSfi}' 



where is given by 



^b.Pb^R^ 



(11) 



(12) 



with Pb being the stationary state of the bath, while the 
system states Puit) are defined by 



~Pi,{t) = TTBi^RMmn)- 



(13) 



We have introduced a set of projectors Hr ~ 
which provides an orthogonal decompo- 
sition of the unit operator [Ib] in the Hilbert space of 
the bath, Y^r^R = ^b, with HrIIr' = UrSr.r'. The 
full set of states jcfl) corresponds to the base where Pb 
is diagonal, which implies J2r^R — Pb- 

It is easy to realize that = V. In physical terms, 
this projector takes in account that each bath-subspace 
associated to the projectors induces a different sys- 
tem dynamics, each one represented by the states Pi^{t). 
Each sub-space can be seen as a sub-reservoir. On the 
other hand, notice that the standard projector Vpx{t) = 
TiBiPri^)} ® Pb = Psi^) ® Pb Ul]; recuperated when 
all the states PR{t) have the same dynamics. Therefore, 
it is evident that the definition of the projector Eq. (fTT|) 
implies the introduction of a generalized Born approxi- 
mation , where instead of a uncorrelated form for the 
total system-environment density matrix, it is assumed a 
classical correlated state. 



By using that ^r = Ib, the system density matrix 
can be written as 

Psit) = Y,TrB{URMt)IlR}^^^ (14a) 
= TrB{PpT{t)} = Y.~PR(*) (14b) 

R 

This equation defines the system state as a sum over 
the states PR{t). Notice that the second line follows from 
the definition of the objects that define the projector 

Eq. In]). 

By writing the evolution Eq. pop in an interaction rep- 
resentation with respect to Hs + Hb, and splitting the 
full dynamics in the contributions Vpj'{t) and Qprp{t), 
where Q = 1 — up to second order in the interaction 
Hamiltonian it follows llli 



Jt 



= f dt'VCT{t)CT{t')VpT{t'), (15) 

Jo 



where Cxit) is the total Liouville superoperator in a in- 
teraction representation. For writing the previous equa- 
tion, we have assumed Qp^iO) = 0, which implies the 
absence of any initial correlation between the system and 
the bath, Pt(S^) = Psi^) ® Pb- Then, the initial condition 
of each state PR (t) can be written as 



PRio) = PrMo). 



(16) 



The parameters Pr are defined by the weight of each 
sub-reservoir in the full stationary bath state 

Pr = TiBiER} = TtbITIrPb} - J2 (^rIPbI^r), (17) 

which trivially satisfies 'YrPr. — 1- 

Now, we split the interaction Hamiltonian as 



(18) 



R,R' 



R,R' 



We notice that when HrHiUri — for R ^ i?', the inter- 
action Hamiltonian can be written as a direct sum Hi = 
Hi^ ®Hi^---® Hi^ © F/^^, • • • , with = n^iJ/HR. 
This case recover the assumptions made in Ref. j22| . 
In fact, without considering the non-diagonal terms in 
Eq. ([5]) [a^'^/ = 0], after a trivial change of notation 
~Pr{^) ~^ PrPr{^) ™ O; llis dynamics reduce to a 
random Lindblad equation. 

In order to proceed with the present derivation, we 
introduce the superoperator identity [U 

[a, [b, •]] = ^[[a, h*] + \{{a, b}+, •}+-(a.S+S.a), (19) 

valid for arbitrary operators d and b. By using this iden- 
tity and the splitting Eq. into Eq. after a 
straightforward calculation the evolution of PR{t) in the 
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Schrodinger representation can be written as in Eq. ([5]). 
The effective Hamiltonians read 

H'Jf = Hs-i- drTTBAlHi, i?/(-r)]ps J. (20) 

The non-diagonal operators De.'r read 

(21) 

while the corresponding superoperators Frr' can be writ- 
ten as 

POO 

Frr'I']^ / dTTi:BAHi^^,{-T)[.]<S>PB„Hi^,„+h.c.). 
Jo 

(22) 

The diagonal contributions follows from the previous ex- 
pressions as Dr = Drr, and Fr[»] = Frr[»]. Further- 
more, we have defined Tr^^l*} = Ttb{^r • H^} and 



(23) 



Notice that these objects correspond to the stationary 
state of each sub-reservoir. 

In obtaining Eqs. to ([^^ we have introduced a 
standard Markovian approximation P, which allows 
to obtain local in time evolutions for the set {pR{t)}, as 
well as to extend the time integrals to infinite. This ap- 
proximation applies when the diagonal and non-diagonal 
correlations of the different sub-reservoirs define the 
small time scale of the problem. In order to clarify the 
introduction of the Markov approximation, we assume 
that the interaction Hamiltonian can be written as 



Hi = 



(24) 



where the operators Va and Ba act on the system and 
bath Hilbert spaces respectively. By using Hj = Hj, the 
previous expressions Eqs. (I^T]) and read 



1 poo 

Dr'H = ^Y1 / dr{x''J^{^r)V^Vf,{~r) + h.c.}, (25) 



and 



poo 

Frr>[>] = / {dTX%%{-T)Vp{-r)[.M + h.c.}. 

(26) 

Here, we have defined the ^^projected bath correlations" 

xfR'i-r) ^ TTB„,{PB^BiTlRB^i~r)}. (27) 

Without taking in account the indexes R and R', this ex- 
pression reduces to the standard definition of bath cor- 
relation P, [3, H, [s^l- Here, the same structure arises 
with projected elements. As the integrals that appears 
in Eqs. (|25p and ([^5]) have the same structure that in the 
standard Born-Markov approximation (3^ . the meaning 
of the previous calculation steps becomes clear. 



Finally, in order to obtain the explicit expressions for 
the matrixes a^'^/ and a'^ , we define a matrix C^^(— r) 
from 



-ffs^^g+-//s = J2 Cfi-,{-T)V^. (28) 



By introducing these coefficients in Eqs. ((25|) and (|26p . 
it is possible to write the operators Dr'r and Frr' [•] as 
in Eq. ([S]). The matrix a^'Jj, is defined by 

pOO 

+ E / dT{xfR>)\-r) C*p,{-t), (29) 

while the diagonal matrix elements follows as a B^ — ^ RR' 
Consistently, without taking in account the indexes R 
and i?', this matrix structure reduce to that of the stan- 
dard Born-Markov approximation [s^ ] . 



Quantum master equation for a system influencing its 
environment 

In Ref . fToj , Esposito and Gaspard deduced a quantum 
master equation intended to describe physical situations 
where the density of states of a reservoir is affected by the 
changes of energy of an open system. While this physical 
motivation is different to that of the GBMA (or in 
general, to the correlated projector techniques \XA), here 
we show that both formalisms can be deduced by using 
the same calculations steps. Therefore, the evolution of 
Ref. 'T^l can also be written as a Lindblad rate equation. 

In Ref. [l3|, the system evolution is derived by taking 
in account the effect of the energy exchanges between 
the system and the environment and the conservation of 
energy by the total (closed) system-reservoir dynamics. 
These conditions are preserved by tracing-out the bath 
coherences and maintaining all the information with re- 
spect to the bath populations. Therefore, the system 
density matrix is written in terms of an auxiliary state 
that depends parametrically on the energy of the envi- 
ronment, which is assumed in a microcanonical state. By 
noting that in the GBMA there not exist any coherence 
between the different sub-reservoirs [see Eg . ([TT]) ]. we re- 
alize that the dynamics obtained in Ref. [l^] can be recov- 
ered with the previous results by associating the discrete 
index R with a continuos parameter e, which label the 
eigenvalues of the reservoir, joint with the replacements 



R 



de n{e), 



(30) 



where n{e) is the spectral density function of the reser- 
voir. Consistently, the system state [Eq. ([H]) ] is written 
as 

Ps{t)^ / de n{e)p{e;t) = / de p{e;t). (31) 
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As in the GBMA, the evolution of p(e; t) can be written 
as a Lindblad rate equation defined in terms of the matrix 
structure Eq. dm with the replacement X^/{/(— t) 
Xee'i-^)' where 

x:^,{^r) = (e'l Bi \e) {e\ Bp \e') exp[-^(e - e')T]. (32) 

This last definition follows from the microcanonical state 
of the reservoir [p^ — )■ 1]. Finally, by introducing the 
matrix elements 

PsAe;t) = {s\-p{e-t)\s'), (33) 

where {|s)} are the eigenstates of the system Hamilto- 
nian, Hs\s) = es\s), the master equation of Ref. fToj is 
explicitly recovered. Due to the energy preservation con- 
dition, in general the evolution involves a continuos para- 
metric coupling between the matrix elements Pss'{£]t) 
and Pss'{£ ± ^;i), where A is a energy scale that char- 
acterize the natural transition frequencies of the system 

We remark that the difference between both ap- 
proaches relies on the assumed properties of the environ- 
ment. In the context of the GBMA, the index R label 
a set of Hilbert subspaces each one defined in terms of 
a manifold of bath eigenstates able to induce, by itself, 
a Markovian system dynamics. Therefore, by hypothe- 
sis, the complete environment does not feels the effects of 
the system energy changes. On the other hand, the ap- 
proach of Esposito and Gaspard applies to the opposite 
situation where, by hypothesis, the density of states of 
the environment vary on a scale comparable to the sys- 
tem energy transitions. The stretched similarity between 
both approaches follows from the absence of coherences 
between the different (discrete or continuous) bath sub- 
spaces. In both cases the system evolution can be written 
as a Lindblad rate equation. 

B. Composite environments 

The previous analysis relies in a bipartite system- 
environment interaction described in a GBMA. Here, we 
arrive to a Lindblad rate equation by considering com- 
posite environments, where extra degrees of freedom U 
modulate the interaction (the entanglement) between a 
system S and a Markovian reservoir B pll | . This formu- 
lation allows to find the conditions under which Eq. ([5]) 
defines a completely positive evolution. 

The total Hamiltonian reads 

Ht = Hs + Hu+Hsu + Hb+Hi. (34) 

As before, Hs represent the system Hamiltonian. Here, 
Hb is the Hamiltonian of the Markovian environment. 
On the other hand, Hjj is the Hamiltonian of the 
extra degrees of freedom that modulate the system- 
environment interaction. The interaction Hamiltonian 
Hi couples the three involved parts. We also consider 



the possibility of a direct interaction between S and U, 
denoted by Hsu- 

As i? is a Markovian reservoir, we can trace out its 
degrees of freedom in a standard way Therefore, 
we assume the completely positive Lindblad evolution 

= :^[Hc,Pcit)] ~ {Dc,Pcit)}+ + Fc[Pc{t)l 

(35) 

with the definitions 
Dc^lY.^^,A]A, Fc[»]^J2^^,A*A]. (36) 

The matrix P(j{t) corresponds to the state of the "com- 
pose system" S-U with Hilbert space He = Hs ® Hu- 
The sum indexes i and j run from one to 1 to (dimTYc)^, 
with dim Tic = diniHs dmiHu- Consistently, the set 
{Ai} is a base of operators in Hci and bij is an arbi- 
trary Hermitian semipositive matrix. 

In order to get the system state it is also necessary 
to trace out the degrees of freedom U. In fact, Psit) = 
Try{pp(t)}, which deliver 

Psit) = Tru{pcit)}^J2(^\Pcii)\R)^ 

R 

R. 

where {|-R)} is a base of vector states in Hu- We no- 
tice that here, the sum structure Eq. (jj]) have a trivial 
interpretation in terms of a trace operation. 

By assuming an uncorrelated initial condition P(j{Q) — 
Ps{0)®Pu(Q), where ^^(0) and /3^(0) are arbitrary initial 
states for the systems S and U, from Eq. ([57]) it follows 
the initial conditions P/^XO) = P^psiO), where 

Pr = {R\Pu{0)\R). (38) 

Therefore, here the weights Pr corresponding to Eq. ([6]) 
are defined by the diagonal matrix elements of the initial 
state of the system U. From now on, we will assume that 
the set of states correspond to the eigenvectors 

basis of Hu, i.e., 

Hu \R) - SR \R) . (39) 

The evolution of the states PR{t) = {R\p(j{t) \R) can 
be obtained from Eq. (|35p after tracing over system U. 
Under special symmetry conditions, the resulting evolu- 
tion can be cast in the form of a Lindblad rate equation, 
Eq. (O. In fact, in a general case, there will be extra con- 
tributions proportional to the components {R\ P(j{t) \R') . 
By noting that 

TTs[{R\pc{t)\R')]^{R\Puit)\R'), (40) 

where pjj(t) — Tisipcit)} is the density matrix of the de- 
grees of freedom U, we realize that the evolution of pR{t) 
can be written as a Lindblad rate equation only when the 
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evolution of pu (t) does not involves coupling between the 
populations {R\ Pu{t) \R) and coherences {R\Pu(i) \R') , 
R 7^ i?', of system U. As is well known [U, M this 
property is satisfied when the dissipative evolution of 
Pij(t) can be written in terms of the eigenoperators 
of the unitary dynamic, i.e., [Hu,Lu\ — louLu- In what 
follows, we show explicitly that this property is sufficient 
to obtain a Lindblad rate equation for the set of matrixes 

{pRit)}- 

First, we notice that the Hamiltonian Hq in Eq. ([55]) 
must to have the structure 

Hc^Hs + Hu+J2Va®L^, (41) 

a 

where Lq are the eigenoperators with a null eigenvalue, 
i.e., [HjjjLq] — 0. With this structure, the populations 
and coherences corresponding to U do not couple be- 
tween them. Therefore, the effective Hamiltonian H'^^^ va 
Eq. (O reads 

H'^^f ^Hs + Y.{R\L^\R)V^. (42) 

a 

After taking the operator base in Tic — Ti-s ® Hu as 
{A,}^{V^^Lu}, (43) 
the superoperators Eq. ([M]) can be written as 

Dc = ^Y.^u2V^^lV^^u, (44a) 

Of, 7 
u,v 

Fc['] = ^C^l/aL^.l/^Lt. (44b) 

a, 7 

With these definitions, by taking the trace operation over 
the system U in the evolution Eq. (I35p . we notice that 
the evolution of the set {pj^it)} can be cast in the form 
of a Lindblad rate equation if the conditions 

J2 {R"\ Ll \R) {R\ i„ \R') = 5r,^r„ all, (45) 

are satisfied. The factor Sb/^rh guarantees that the 
evolution of the set {pjiit)} do not involve the terms 
{R\ Pcit) \R') , R 7^ R', and in turn implies that the pop- 
ulations and coherences of U do not couple between them. 
On the other hand, defines the matrix elements cor- 
responding to the structure Eq. ([S]). The diagonal con- 
tributions follows from Eq. (|45p by taking R — R'. 

The set of conditions Eq. (|45p can be simplified by 
taking the base 

Lu ^ \n'){nl (46) 

which from Eq. ((39| satisfy [Hu^Lu] = {s-jz — e-R,')-L„. 
Thus, Eq. (|45p can be consistently satisfied if we impose 

bZ2 = 0, for u^v. (47) 



After changing J2u^J2n,TZ' i'^ 63; we get 

^BM' — "{RM'){RM')^ '^R " °(_R,i?.)(_R,fl) ' \^^) 

where we have used that TZ and TZ' are dumb indexes. 
This result demonstrate that the evolution induced by 
the composite environment can in fact be written as a 
Lindblad rate evolution Eq. ([5]) with the matrix elements 
defined by Eq. (p|. 

From our previous considerations we deduce that Lind- 
blad rate equation arise from microscopic tripartite inter- 
actions having the structure 

Hi = La® HsB + Y,U®H^B + LI ® (i?5s)^ (49) 

U 

where [Hu,Lo] = 0, and L„ ^ \R){R'\ with R ^ R' . 
On the other hand, Hg^ are arbitrary interaction terms 
between the system 5* and the Markovian environment 
B. In fact, the structure Eq. (|49)) guarantees that the 
populations and coherences of U do not couple between 
them, which in turn implies that the evolutions of the 
system S is given by a Lindblad rate equation. 

Completely positive condition 

We have presented two different microscopic interac- 
tions that lead to a Lindblad rate equation. In order 
to use these equations as a valid tool for modeling open 
quantum system dynamics it is necessary to establish the 
conditions under which the solution map ^5(0) — > Ps{i) 
is a completely positive one. For an arbitrary Lindblad 
rate equation this condition must to be defined in terms 
of the matrixes a^'Jj/ and . 

In order to find the allowed matrix structures, we no- 
tice that the evolution Eq. (|35p is a completely positive 
one when hij ^Ir r')(r w) ^ semipositive defined 
matrix. Therefore, by using Eq. (|48p we arrive to the 
conditions 

Kr\,\>0. I4''I>0, Vi?,i?', (50) 

i.e., for any value of R and R' both kind of matrixes must 
to be semipositive defined in the system indexes a, 7. The 
condition \a'^\ > has a trivial interpretation. In fact, 
when a'^ff, — 0, there not exist any dynamical coupling 
between the states Pji{t). Thus, their evolutions are de- 
fined by a Lindblad structure that under the constraint 
> define a completely positive evolution. 

C. Quantum random walk 

By using the similarity of Eq. ([5]) with a classical rate 
equation [33[ , here we present a third derivation by con- 
structing a stochastic dynamics that develops in the sys- 
tem Hilbert space and whose average evolution is given 
by a Lindblad rate equation. 
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First, we assume that the system is endowed with a 
classical internal degree of freedom characterized by a set 
{R} of possible states. The corresponding populations 
PR,{t) obey the classical evolution 



1^ iH'nPnit) 

R' 
R'^R 



E ^RR'PR'ii)^ (51) 

R' 
R'^R 



with initial conditions Pr{Q) = Pr, and where the co- 
efficients {'^fjifj} define the hopping rates between the 
different classical states R. 

To each state R we associate a different Markovian 
system dynamics, whose evolution is generated by the 
superoperator 



Cb=C 



H 



(52) 



with -C_ff[»] = (^*/^)[-ffsj •] and a standard Lindblad 
contribution [•] = —{Dii,»}++F]i[»]. Therefore, each 
state R defines a propagation channel with a different 
self-dynamic. The system state follows by tracing out 
any information about the internal state. Thus, we write 



(53) 



where each state Pji{t) defines the system state given that 
the internal degree of freedom is in the state R. Consis- 
tently, the initial condition of the auxiliary states reads 

~Pr{^) = PrPs{^)- 

Finally, we assume that in each transition i? — > i?' of 
the internal degree of freedorn, it is applied a completely 
positive superoperator En 0, H, 0] , which produces a dis- 
ruptive transformation in the system state. 

The stochastic dynamics is completely defined after 
providing the self-channel dynamics, defined by {Cji}, 
the set of rates {7^^/^} and the superoperators {£r}- By 
construction this dynamics is completely positive. The 
explicit construction of the corresponding stochastic re- 
alizations, which develop in the system Hilbert space, is 
as follows. When the system is effectively in channel i?, 
it is transferred to channel R' with rate ^pjn- Therefore, 
the probability of staying in channel R during a sojourn 
interval t is given by 



PfHO=expH 



Ir'rV 



(54) 



R' 
R'^R 



This function completely defines the statistics of the time 
intervals between the successive disruptive events. As in 
standard classical rate equations, when the system "jump 
outside" of channel R, each subsequent channel R' is se- 
lected with probability 



tR'R. — 



iR'R 



R" 1R"R 
R"^R 



(55) 



in such a way that X]_r' ^^i'^? = 1- Furthermore, each 
transference R ^ R' , is attended by the application 



of the superoperator f^, which produces the disrup- 
tive transformation Pfi(t) — » £fi[pj^{t)]. This transformed 
state is the subsequent initial condition for channel R'. 

The average over realizations of the previous quantum 
stochastic process, for each state Pfj{t), reads 



PR{t) = Prit)e'^"p{0)+ drPl,^\t-r)e^'-^->^- 



lRR'^R'[PR'i^)h 



(56) 



R' 
R'^R 



The structure of this equation has a clear interpretation. 
The first contribution represents the realization where 
the system remains in channel R without happening any 
scattering event. Clearly this term must be weighted by 
the probability of not having any event in the time inter- 
val {t, 0), i.e., with the probability PQ^'\t). On the other 
hand, the terms inside the integral correspond to the rest 
of the realizations. They take in account the contribu- 
tions that come from any other channel i?', arriving at 
time r and surviving up to time t in channel R. Dur- 
ing this interval it is applied the self-channel propagator 
exp[(i — t)jCii]. As before, this evolution is weighted by 
the survival probability PQ^\t — r). 

By working Eq. (|56p in the Laplace domain, after a 
simple calculation, it is possible to arrive to the evolution 

J^pRit) = Y^Hs,PR{t)] - {DRrPR{t)}++FR[pj,{t)] 

- E lR'RPRit)+ E iRR'^R'lPR'it)]- 
R' R' 
R'^R R'^R 

(57) 

We notice that this expression does not corresponds to 
the more general structure of a Lindblad rate equation, 
Eq. ([5]). Nevertheless, there exist different non-trivial 
situations that fall in this category. As we demonstrate 
in the next section, the advantage of this formulation is 
that it provides a simple framework for understanding 
some non-usual characteristics of the system dynamics. 



III. NON-MARKOVIAN DYNAMICS 

In this section we obtain the master equation that de- 
fine the evolution of the system state Ps{t) associated to 
an arbitrary Lindblad rate equation, Eq. ([S]). 

In order to simplify the notation, we define a column 
vector defined in the i?-space and whose elements are the 
states pj^, i.e., \p) = {p^, P2, ■ ■ ■ Pr, ■ ■ ■)'^ , where T denote 
a transposition operation. Then, the evolution Eq. ([5]) 
can be written as 



dm) 

dt 



CH\pit))+M\pit)). 



(58) 
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where 

of M reads 



-{i/K)[Hs,»], and the matrix elements 



^B.R' [•] = 5r,r, { -J-[H'r, •] - {Dr, .}+ + Fr,[»] 



'Frrj [.] - Sr^rj {Dr„r, .} + , (59) 



R" 



where H'j^ = H'j^J — Hs, is the shift Hamiftonian pro- 
duced by the interaction with the reservoir. The initial 
condition reads \p{0)) = \P) Ps(0), where we have intro- 
duced the vector \P) — (Pi, P2, ■ ■ ■ Pr, ...)"'". The system 
state Eq. ^ reads Ps{t) = (1 \p{t)) , where |1) is the 
row vector with elements equal to one. Notice that due 
to the normalization of the statistical weights it follows 
(1|P) = 1. 

From Eq. the system state can be trivially written 
in the Laplace domain as 



Ps(w) = (1 



u- (Ch+M) 
(l|G(u)|P)ps(0) 



l^)Ps(O), 



(60a) 
(60b) 



where u is the conjugate variable. Multiplying the 
right term by the identity operator written in the form 
1/ (1| G(m)[u - {jCh + M)] \P) , it is straightforward to 
arrive to the non-local evolution 

= CnlPsit)] + £ drUt ~ r)[psir)], (61) 

where the superoperator L(t) is defined by the relation 



(1|(G(m)M|P) 



(l|(G(w) |P)L(w)[. 



(62) 



In general, depending on the underlying structure, the 
evolution Eq. (j6ip involves many different memory ker- 
nels, each one associated to a Lindblad contribution. 

We notice that a similar master equation was obtained 
in Refs. [2l|, [23|. Nevertheless, here the dynamics may 
strongly departs with respect to the evolutions that arise 
from Lindblad equations with a random rate [a^}j, — 0] . 
In fact, the previous calculation steps are valid only if 



lim„_^o {l\u<G{u) \P) = 0. 



(63) 



By using that limf^oo /(i) — ^^''^-u^o uf{u), this con- 
dition is equivalent to limt^oo (1| G(i) |P) — 0. In the 
general case a^'Jj/ 7^ 0, Eq. ([55)1 is not always satisfied. 
In this situation, the density matrix evolution becomes 
non-homogenous and the stationary state may depends 
on the system initial condition. In general, this case 
may arises when the diagonal contributions are null, i.e., 
"fl^ = and a^"^, ^ 0. We remark that these matrix 
structures values are completely consistent with the con- 
ditions Eq. (|5D|) . On the other hand, in the context of 
the GBMA, this case arise when the diagonal sub-bath 



correlations are null, x'tfR{~''') = Oj which in turn im- 
plies that the interaction Hamiltonian Eq. ^TE\i satisfies 
UrHiUr, = if i? = i?'. 

In order to characterize the dynamics when the condi- 
tion Eq. is not satisfied, we introduce the difference 



Spsiu) 



= psiu) ~ - lim (l|u(G(u) |P)/95(0),(64a) 
u ti^o 

= (11 <G(u) - - lim u<G(u) IP) ps(0), (64b) 

U u^O 



{l\dGiu)\P)psiO), 



(64c) 



where now the pseudo-propagator 6G{u) satisfies 
lim 1,^0 (l|u(5'G(u) |P) = 0. Therefore, Sps{t) satisfies an 
evolution like Eq. ([CT]) where the kernel is defined by 
Eq. dMl) with G(u) 6G,{u). Notice that the system 
state, even in the stationary regime, involves the contri- 
bution lim„^o (1| m(G(m) |P) pg(0), that in fact depends 
on the system initial condition. 

In the next examples we show the meaning of this prop- 
erty, as well as its interpretation in the context of the 
stochastic approach. 



A. Dephasing environment 

Here we analyze the case of a qubit system interacting 
with a dispersive reservoir [J, [l5| whose action can be 
written in terms of a dispersive Lindblad rate equation. 
We assume a complex reservoir with only two subspaces, 
R — a,b, whose statistical weights [Eq. (fT7| ] satisfy Pa + 
Pfc = 1. Thus, the system state reads 



(65) 



A generalization to an arbitrary number of sub-reservoir 
is straightforward. 

The evolution of the auxiliary states are taken as 



dt 

d 
dt 



Pait) = -laiPait) -'^zPai'tyz] 

-IbaPaii) + Jab<^zPb{t)'^z, (66a) 

Pbit) = -IblPbit) - ^zPbit)crz] 

-labPbit) +lba^zPa{i)(^z, (66b) 



where cr^ is the z Pauli matrix. The completely positive 
conditions Eq. ([50| imply 



(67a) 
(67b) 



la > 0, 7fc>0, 
lab > 0, 7fca > 0. 

By denoting the matrix elements by (P = a, 6) 

. Kit) $+(0 \ 



(68) 
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the evolution corresponding to the populations read 
d 



^n^(t) = -76an^(i)+7a6nf(t), (69a) 



^n±(t) = -7„,n±(i) + 7,„nJ(t), (69b) 

with n^(0) = PRn|(0), while for the coher ences we ob- 
tain 



<{t) = -{la+lba)^a{t)--lab^i{t), (70a) 



d 

dt 

J^'^bit) = Hlb + lab)^bit)-lba^tit), (70b) 

with $±(0) = Pfl«>|(0). For expressing the initial condi- 
tions we have trivially extended the notation Eq. ([55]) to 
the matrix elements of Psit)- 

We notice that all coherences and populations evolve 
independently each of the others. From the evolution of 
the populations Eq. ([69|) it follow 



FIG. 1: Normalized coherences $s(t)/*s(0) = h{t), Eq. ffi)! . 
In the upper curve the parameters are — 0.1, 7;, = 1, and 
lab = 76a = 0. In the lower curve they are 7^ = 0.1, 7^ = 1, 
7„(, — 1, and 75^ = 0.1. The rates are expressed in arbitrary 
units (a.u.). In both curves we take Pa = 0.1 and Pi, = 0.9. 



dt 
d_ 
dt 



Tr[p„(t)] = -7fcaTr[p,(t)]+7ahTr[p,(i)], (71a) 
Tr[p,(i)] = -7afcTr[p,(t)]+76aTr[p,(t)], (71b) 



with Tr[p„(0)] + Tr[pf,(0)] = P„ + Pb = 1, which implies 
that the trace of the auxiliary states perform a classical 
random walk. 

From Eqs. (|65p and ((69l) it becomes evident that the 
populations of the system remain unchanged during all 
the evolution. On the other hand, the dynamic of 
the coherences can be obtained straightforwardly in the 
Laplace domain. From Eq. ((TOl) we get 

^t{u) = /la6(u)$|(0), <i>±(w) = /ib,(w)<i>|(0), (72) 

where we have introduced the auxiliary function 

^ iPa-Pbhab + Pa{u + -/b) 

Ibaiu + 7a) + labiu + 7b) + (" + 7a)(w + 76) ' 

(73) 

Therefore, from Eq. (|65p the matrix elements of pg(t) 
read 

n|(0 - n±(o), $!(<) - /i(i)<i>t(0), (74) 

where h{t) = hab{t) + hba{t), gives the coherences decay. 
From these solutions, it is straightforward to obtain the 
corresponding system evolution 



dp sit) 



dt 



dTKit^T)C[ps{T)], 



(75) 



with £[•] = (— • +az »<Jz) and K{u) = [1 — uh{u)]/h{u). 

In order to check the completely positive condition, we 
write the solution map as 



with g±{t) = [1 ± h{t)]/2. This maj 
positive at all times if g±{t) > 
implies the constraint 

\h{t)\ < 1. 



lappin 



is completely 
1, and in turn 



(77) 



In the upper curve of Fig. ([T]) we plot the normalized 
coherences $g(i)/$g(0) — h{t) for the case in which 
the non-diagonal rates are null, 7^^ = ^^^^ — 0. Then, 
the dynamics reduce to a superposition of exponential 
decays, each one participating with weights Pa and P^. 

In the lower curve of Fig. ([T|) the non-diagonal rates 
are non-null, while the rest of the parameters remain the 
same as in the upper curve. In contrast to the previous 
case, here the coherence decay develops an oscillatory 
behavior that attain negative values. Clearly, this regime 
is unreachable by a superposition of exponential decays. 

In both cases, the condition Eq. (1771) is satisfied, guar- 
anteeing the physical validity of the respective solutions. 



Stochastic representation 



Psit) ^ g+{t)p{0) + g-{t)azp{0)az 



(76) 



The evolution Eq. (|66)) admits a stochastic interpreta- 
tion like that proposed previously. The stochastic tra- 
jectories can be simulated with the following algorithms. 
First, for being consistent with the initial condition, the 
system initialization must be realized as follows 

i) Generate a random number r S (0, 1). 

ii) li r < Pa {r > Pa) the dynamic initialize in channel 
a (6) with pM = ps{Q) [pM = psm- 

Trivially, with this procedure the channel a ih) is ini- 
tialized with probability Pa (Pb)- 

By comparing Eqs. (|66p and (|57p . the scattering super- 
operator results £[•] = az • a^, which does not depends 
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on the channel (a and b) . It action over an arbitrary state 
[Eq. dSS])] is {R = a, b) 



Kit) 



(78) 



Therefore, its application implies a change of sign for 
the coherence components. On the other hand, the self- 
dynamics Eq. ([5^ of each channel is defined by Ca/b [•] = 

7a/b(-»+0-^»cr2)- 

With the previous information, the single trajectories 
can be constructed with the following algorithm: 

1) Given that the system has arrived at time ti to chan- 
nel a, generate a random number r £ (0, 1) and solve for 

(^i+i — ti) from the equation Pg°''(tj+i — ti) ~ r, where 

Po'"'(i) = exp[-7,,t]. 

2) For times satisfying t e (ii+ijti), the dynamics 
in channel a is defined by its self- propagator, p^(i) — 

eyip[{t-t,)Ca]Pait^). 

3) At time t^+i the system is transferred from channel a 
to &, implying the transformation pj,(ti+i) — > £[p^{ti^i)] 
and the posterior resetting of channel a, defined by 
p„(i.+i)^0. 

4) Go to 1) with a b and i ^ i + 1. 

At this point, it is immediate to realize that the clas- 
sical rate equations Eqs. ([55]) and (|7ip arise straightfor- 
wardly from the (transfer) jumps between both channels. 
The corresponding stationary traces read 



Tr[Pa(c 



lab 



lab + Iba 



Tr[Pb(oo)] 



Iba 



lab + Iba 



(79) 



which do not depend on the system initial state. 

In contrast with the population evolution, some non- 
standard dynamical properties can be found in the coher- 
ences evolution when la ~ lb — 0- Fig. ^ we show 
the normalized coherences ^s(t)/^^(0) — h{t) corre- 
sponding to this case. In the inset, it is shown a typical 
stochastic realization of the coherences of the auxiliary 
matrixes p^it) and Pb{t) obtained with the previous al- 
gorithm. As expected, in each application of £ the co- 
herences are transferred between both channels with a 
change of sign. We also show an average over 500 real- 
izations. We checked that by increasing the number of 
realizations, the average behavior result indistinguishable 
with the dynamics Eq. ((74l) . 

In strong contrast with the previous figure, in Fig. Q 
the stationary values of the coherences are "not null and 
depend on the initial condition." In fact, their normal- 
ized asymptotic value is limt_^oo ^s(0/^s(0) — —0.654. 
This characteristic is consistent with the breakdown of 
condition Eq. (|63p and can be understood in terms of 
our previous analysis. By taking = 7^ = in Eq. (I72p 
we get 



Pa{u + lab)-Pblab;u± 



U[U + lab + lba\ 



*s(0), 



(80) 



FIG. 2: Normalized coherences $|(t)/$|(0) = h{t), Eq. ffi)! . 
The parameters are 7a = 76 = 0, 7a6 = 1. Iba = 0.1, with 
the statistical weights Pa = 0.1 and Ph = 0.9. The noisy 
curve correspond to an average over 500 realizations of the 
trajectories defined in the text. The inset show a particular 
realization for the coherences $ J (t) and (t) of the auxiliary 
matrixes p^{t) and Pi,{t) respectively. 



which implies the asymptotic value 



hm $±(t) 

E— i-OO 



[Pa - Pb 



lab 



lab + Iba 



*5(0), 



(P,-n)Tr[p„(^)]$|(0). 



(81a) 
(81b) 



This last expression can be easily interpreted in terms 
of the realizations of the proposed stochastic dynamics. 
From the inset of Fig. (I2|), it is clear that, in spite of a 
change of sign, the coherence transferred between both 
channels does not change along all the evolution. In fact, 
notice that due to the election 7^ = 7^ = 0, the self- 
propagators of both channels [see previous step 2)] are the 
identity operator. Therefore, all realizations that begin 
in channel a [measured by Pa] that are found in channel 
a in the stationary regime (measured by Tr[p^(oo)]), con- 
tributes to the stationary value of the coherence <I>^ (t) 
with the value $^(0). This argument explain the con- 
tribution proportional to Pa Tr[p„(oo)]$|(0) in Eq. ((5T|) . 
On the other hand, a similar contribution is expected 
from the realizations that begin in channel b. Neverthe- 
less, due to the action of the superoperator £ [Eq. (|78p ] 
they contributes with the opposite sign. 

By adding the contributions of both auxiliary matrixes, 
from Eq. ([5T|) the stationary system coherences reads 



lim $±(i) 



(Pa - Pb, 



Tab Tba 
lab + Iba 



$1(0)^0, (82) 



This expression fits the stationary value of Fig. ([2]). 

The stochastic realizations corresponding to the sys- 
tem coherence ^g{t) can be trivially obtained from the 
the realizations of $0 (i) and ^^(t). By adding the up- 
per and lower realizations of the inset of Fig. ^ , we get 
a function that fiuctuates between the values ±$|^(0). 
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By considering the initial conditions and the superoper- 
ator action from these reahzations it is also possible to 
understand the four contribution terms of Eq. ((82|) . Fi- 
nally, we remark that when any of both channels have 
a non-trivial self-dynamics, the coherences vanish in the 
stationary regime, losing any dependence on the system 
initial condition ^^(O) [see Fig. ([T])]. 



B. Depolarizing reservoir 

Another example that admits a stochastic representa- 
tion is the case of a depolarizing reservoir [J, [l^ , which 
is defined by the superoperator 



£[•] = {cTx • CT^E + CTy • cry)/2 



(83) 



where and ay are the x and y Pauli matrixes respec- 
tively. For simplifying the analysis we assume channels 
without self-dynamics. Therefore, the evolution reads 



where nj]~(c») and II~{qo) follows after changing a <-> 6. 
This result has an immediate interpretation in the con- 
text of the stochastic approach. In fact, the last frac- 
tional factors correspond to the "natural" stationary so- 
lutions of Eq. ([86)1 . This solution is corrected by the 
terms in brackets, which in fact take in account the 
system initialization [notice that n+(0) + n^(0) ^ 1] 
and the transformations induced by the superoperator 
£ Eq. l|85|) . Finally, the system stationary populations 
nf (oo) = n±(oo) + n^(oo) reads 

n| (OO) = n| (0) P-^-b + Pblba ^ (O) P-^ba+Pblab 
^ ^ lab+lta ^ lab+lba 

(89) 

As in the previous case, the dependence of the stationary 
state in the initial conditions is lost when the channels 
have a proper dissipative self-dynamics. 



IV. SUMMARY AND CONCLUSIONS 



J^Pait) - ~lba~Pait)+lab£[~Pb{t)], (84a) 
J^Pbit) = -labPbit)+lba£[Pam- (84b) 

The action of the superoperator £ over the states Puit) 
[Eq. is given by (R = a, b) 



£[PR{t)\-y n+(t) 



(85) 



Therefore, its application destroy the coherences compo- 
nents and interchange the populations of the upper and 
lower states. 

The populations of the auxiliary states evolve as 



di 
d_ 

di 



Kit) 



-lbaKit)+labK(t), 



^;it) = -labKit)+lbaKit). (86b) 



subject to the initials conditions n+(0) — Pall^{0) and 
n^(0) = PbUg{0). The evolution of 11+ {t) and n-{t) 
follows after changing a ^ b. Notice that this splitting 
of the population couplings follows from the superoper- 
ator action defined by Eq. ([SS)) . On the other hand, the 
coherences evolution read 



d 
di 



'^tit) = -J,a^iit)^ -<i>tit) = -7.,1>^(i). (87) 



Therefore, in this case the stationary coherences are null. 
This fact also follows trivially from Eq. ([55]) . In contrast, 
the stationary populations reads 



n+(^) = [n+(o)F, + n^(o)p,] 
n-(c») = [n+{Q)Pa + Us{o)Pb] 



lab 



lab + Iba 
Iba 

lab + Iba 



, (88a) 
, (88b) 



We have presented a new class of dynamical master 
equations that provide an alternative framework for the 
characterization of non-Markovian open quantum system 
dynamics. In this approach, the system state is written 
in terms of a set of auxiliary matrixes whose evolutions 
involve Lindblad contributions with coupling between all 
of them, resembling the structure of a classical rate equa- 
tion. 

We have derived the previous structure from different 
approaches. In the context of the GBMA, a complex 
structured reservoir is approximated in terms of a direct 
sum of Markovian sub-reservoirs. Then, the Lindblad 
rate structure arises by considering arbitrary interaction 
Hamiltonians that couple the different subspaces associ- 
ated to each sub-reservoir. The matrix structures that 
define the system evolution are expressed in terms of the 
projected bath correlations. 

On the other hand, we have derived the same structure 
from composite environments, where the entanglement 
between the system and a Markovian environment is 
modulated by extra unobserved degrees of freedom. The 
Lindblad rate structure arises straightforwardly when the 
tripartite interaction Hamiltonian that involve the three 
parts does not couple the coherences and populations of 
the extra degrees of freedom. This scheme also allows 
to find the conditions under which an arbitrary Lindblad 
rate equation provides a completely positive evolution. 

Due to the apparent similarity of the evolution with a 
classical rate equation, we have also formulated a quan- 
tum stochastic dynamics that in average is described by a 
Lindblad rate equation. The stochastic dynamic consists 
in a set of transmission channels, each one endowed with 
a different self-system evolution, and where the transi- 
tions between them are attended by the application of 
a completely positive superoperator. This formalism al- 
lows to understand some amazing properties of the non- 
Markovian dynamics, such as the dependence of the sta- 



tionary state in the initial conditions. This phenomenon 
arise from the interplay between the initial channel oc- 
cupations and the structure of the stochastic dynamics. 
We exemplified our results by analyzing the dynamical 
action of non-trivial complex dephasing and depolarizing 
reservoirs over a single qubit system. 

In conclusion, we have presented a close formalism that 
defines an extra class of non-Markovian quantum pro- 
cesses that may be of help for understanding different 
physical situations where the presence of non-local effects 
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is relevant [13, [H, [H, [H, [H H [M HI . 
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